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SOFTWARE FOR THE GALE TRANSFORM OF FEWNOMIAL 
SYSTEMS AND A DESCARTES RULE FOR FEWNOMIALS 

DANIEL J. BATES, JONATHAN D. HAUENSTEIN, MATTHEW E. NIEMERG, 

AND FRANK SOTTILE 


Abstract. We give a Descartes’-like bound on the number of positive solutions to a 
system of fewnomials that holds when its exponent vectors are not in convex position 
and a sign condition is satisfied. This was discovered while developing algorithms and 
software for computing the Gale transform of a fewnomial system, which is our main 
goal. This software is a component of a package we are developing for Khovanskii-Rolle 
continuation, which is a numerical algorithm to compute the real solutions to a system 
of fewnomials. 


Determining the real solutions to a system of polynomial equations is a challenging 
problem with real-world applications. Typically, few of the complex solutions are real and 
even fewer satisfy meaningful constraints such as positivity. Methods based on symbolic 
computation j21j or numerical homotopy continuation (4[ [22] first find all complex solu¬ 
tions. ft is an active area of research to develop methods that find only the meaningful 
real solutions. Among these are exclusion methods mm, methods based on semidehnite 
programming [15] . the critical point method [T21 DU] , and one proposed by two of us [5J 
that is based on the Khovanskii-Rolle Theorem [15] from the theory of fewnomials. 

A system of polynomial equations with few monomials is a fewnomial system. Such a 
system enjoys bounds on its number of real solutions independent of its number of complex 
solutions PEI EES!- Khovanskii-Rolle continuation is a symbolic-numeric method that 
can find all positive real solutions of a fewnomial system without first finding all complex 
solutions. In fact, it is the first algorithm for computing positive solutions with complexity 
depending primarily on the fewnomial bound on the number of real solutions, rather than 
the number of complex solutions or the ambient dimension. The core algorithm was 
described in [5], along with a proof of concept implementation. 

Khovanskii-Rolle continuation rests upon the fundamental notion of Gale duality for 
polynomial systems , which is a scheme-theoretic isomorphism between complete intersec¬ 
tions of Laurent polynomials in an algebraic torus and complete intersections of master 
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functions in the complement of a hyperplane arrangement [8|. The support of the poly¬ 
nomials annihilates the weights of the master functions, whence the term Gale duality. 

The Khovanskii-Rolle continuation algorithm finds all real solutions to a system of 
master functions. Computing approximations to the real solutions of a fewnomial sys¬ 
tem F requires that F be converted into its dual Gale system G which is solved using 
Khovanskii-Rolle continuation, then the numerical approximations to solutions of G are 
converted back to numerical approximations to solutions of the original fewnomial system 
F. Symbolic choices made in constructing the dual system affect the performance of the 
Khovanskii-Rolle continuation algorithm. 

One such choice is that of a basis for the nullspace of the matrix of exponent vectors of 
the monomials of F. When the nullspace meets the positive orthant, a bound lower than 
the fewnomial bound holds. Choosing a positive basis element results in far fewer paths 
to be followed in the Khovanskii-Rolle continuation algorithm. This happens when the 
exponent vectors of the original polynomial system are not in convex position and a sign 
condition holds on the coefficient matrix. Such sign conditions are also found in recent 
work giving very strong bounds on positive solutions mum and are considered to be 
multivariate versions of Descartes’ rule of signs. 

In addition to this new Descartes’-like bound, we make explicit the algorithms for both 
the creation of G from F and the computation of solutions of F from solutions of G. 
These steps were described existentially in [8]. We have implemented these algorithms 
in a software package galeDuality, which is available at each authors’ website. This 
will be the front end for a software package that we are developing for Khovanskii-Rolle 
continuation. 

Section Q] gives a description of Gale duality and sketches Khovanskii-Rolle continuation, 
illustrating it through an extended example. In Section [2j we present our Descartes’-like 
bound for fewnomial systems and establish a key lemma about the location of points to 
be tracked in the algorithm. Explicit pseudocode for using Gale duality to compute mas¬ 
ter functions (Algorithm 13. lj) and to transform approximations to solutions of the Gale 
system into approximations to solutions of the fewnomial system (Algorithm I3.3|l is pro¬ 
vided in Section [31 along with details for portions of these algorithms. One subroutine of 
Algorithm [3J] requiring more care is covered in Subsection 13.31 Finally, a brief description 
of our software package is provided in Section [4j 

1. Khovanskii-Rolle continuation and the Gale transform 

We give an example of Gale duality for polynomial systems and the Khovanskii-Rolle 
continuation algorithm. For a complete treatment, see M. 

1.1. An example of Gale duality. Consider the system of Laurent polynomials, 

v 2 w 3 — lluvw 3 — 33 uv 2 w + 4 v 2 w + 15 u 2 v + 7 = 0 , 

(1.1) v 2 w 3 + 5uv 2 w — Av 2 w — 3u 2 v + 1 = 0, 

v 2 w 3 — lluvw 3 — 31 uv 2 w + 2 v 2 w + 13 u 2 v + 8 = 0 . 
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This system has 20 complex solutions. Six are real and three are positive, 

(1.194,0.374,1.231), (0.431,0.797,0.972), (0.613,0.788,0.850). 
Solving the equations in (11.111 for the monomials v 2 w 3 , v 2 w, and uvw 3 , gives 


v 2 w 3 


= 1 — u 2 v — uv 2 w , 


( 1 . 2 ) 


Since 


v 2 w = \ — u 2 v + uv 2 w , 

§i 


and 


uvw 


= yj(1 + u 2 v — 3uv 2 w) 


(■ uv 2 w ) 3 • (v 2 w 3 ) = u 3 v 8 w 6 = (u 2 v) ■ ( v 2 w ) 3 • (uvw 3 ) and 


(u 2 vy ■ ( v 2 w 3 ) 3 = u 4 v 8 w 9 = ( uv 2 w ) 2 • ( v 2 w ) • ( uvw 3 y , 

we may substitute the expressions on the right hand sides of (II. 2\i for the monomials v 2 w 3 , 
v 2 w } and uvw 3 in these expressions to obtain the system 

(uv 2 w) 3 ■ (l — u 2 v — uv 2 w) = (u 2 v) ■ (| — u 2 v + uv 2 w ) 3 • (jj(l + u 2 v — 3uv 2 w )) 

(u 2 vy ■ (l — u 2 v — uv 2 w ) 3 = ( uv 2 wY ■ (| — u 2 v + uv 2 w) ■ (yj(1 + u 2 v — 3uv 2 w)y . 

Writing x for u 2 v and y for uv 2 w and solving for 0, these become 

g \= y 3 (l-x-y) - x(| - x + y) 3 (g(l + x - 3y)) = 0, and 

/ := x 2 (l -x-y) 3 - y 2 (\-x + y){^{l + x-3y)) 2 = 0. 

This system also has 20 solutions with six real, outside of the five lines where the degree 
one factors vanish. Figure [T| shows the curves the equations (11.31) define and the five lines. 
Three solutions lie in the pentagon. They correspond to the positive solutions of (11.11) . 


(1.3) 



Figure 1. Curves and lines 

This transformation from the original polynomial system (II.ip to the system (II.3ft is 
called Gale duality for polynomial systems. The second system has three forms which are 
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equivalent in the pentagon. We first rewrite (j 1.3 [i as a system of rational master functions 

— T — 77/^fl — T — 7/^ 

(1 ' 4) vHi-x + y) (jy(l + x - 3y)) 2 = x(\-x + yY (±?(1 + z - 3»)) = 

Then we take logarithms to obtain the last form. 

5 21og(x) + 31og(l—x—y) — 21og(y) — log(|—x+j/) — 2 log (^(l+x—3?/)) = 0 

31og(y) + log(l-x-y) — log(x) - 31og(i-a;+2/) - log (^(l+x-3y)) = 0 


1.2. Geometry of Gale duality. Let R x = R/{0} be the nonzero real numbers. Sup¬ 
pose that F : (R x ) n —y R n is given by Laurent polynomials that have a total of n+£ 
distinct monomials, none of which is a constant. For b G R n , the system F(x) = b of 
Laurent polynomials has an equivalent expression. Let A := {cq,..., a n+ g} C Z n be the 
set of n+£ exponents of the monomials in F, which we call the support of F. Consider 
the map 

<p A : (R x ) n —> (R x ) n+£ c R n+£ 

Affine-linear forms L(z) = A pull back to polynomials with support A, and so F(x) = b 
is the pullback of a codimension n (dimension £) affine linear space F£ of R n+£ defined by 
C(z) = b, where C G R nx ( n +D is the coefficient matrix of F. We define Sp to be the set 
of solutions of F(x) = b. Then, 

(1-7) <p A (S F ) = <p A ((R x ) n )nH. 

Let -0: R f —> H be any affine isomorphism. This is given by a point A = -0(0) G H and 
n+£ linear forms L\,, L n+ c on R^ which span the dual space of R , 

(1.8) ip(y) = (Li(y) + Ai, L 2 (y) + A 2 ,..., L n+l (y ) + X n+e ) G H , for y G R f . 

The forms L 1; ..., L n+i span the kernel of C, C(Li(y ),..., L n+ i(y)) = 0. Dehne 

S G :=^-\ip A {S F ))<Z^ 1 . 

This set (which is isomorphic to Sp as a scheme) satisfies a simple system of equa¬ 
tions in Rh 

Any integer linear relation with coefficients $ G Z among the exponents in A , 

n+£ 

o = y , 

i= 1 

gives a monomial equation that is satished on <p u 4 ((R x ) n ) in (R x ) n+ 0 

n+£ 

(1.9) n s = i - 

i— 1 

This pulls back under 0 to an equation involving a rational master function, 

n+i 

(1-10) U(Li(y) + Xit = 1, 

i= 1 
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which is defined in the complement of the hyperplanes given by the affine forms L, (y) + X,. 

Writing elements of A as column vectors yields an n x (n+£) integer matrix, also 
written A. An element of the integer nullspace of A gives a monomial equation (11.9ft 
vanishing on yj_ 4 ((M x ) n ). Choosing a Z-basis for this nullspace gives a set G of £ equations 
involving master functions (II. 10)1 . and these define the set Sq- Assuming there are finitely 
many solutions to F{x) = b in (C x ) n , these i monomial equations (jl.9(1 that cut out the 
torus <+ 4 ((C x ) n ) in (C x ) n+ ^ restrict to a complete intersection in H and any r of these 
master functions defines a subset in the hyperplane complement of codimension r. We 
call this set G of £ master functions a Gale system dual to the system F(x) = b of Laurent 
polynomials. We summarize some results in [7, <3]. 

Proposition 1.1. Let F(x) = b be a system of n Laurent polynomials on (M x ) n with 
support A consisting of n+£ nonzero exponents of monomials expressed as a matrix 
A G Z, nx( ' n+e ) and finitely many zeroes in (C x ) n . Necessarily the exponent vectors span 
a full rank sublattice ZA of Z n . Let be the map (1 1.6 j) and C G R nx ( n +P be the coeffi¬ 
cient matrix of F. The kernel of is the Pontryagin dual to 7L n /7LA. 

A Gale system dual to F(x) = b is given by a choice of a basis B = {/3^,... ,/3^} C 
Z n+£ for the integer nullspace of A, a choice of a basis L := {Li ,..., L n+ e} for the kernel 
of C T , and a choice of a point A G K n+£ with C T A = b. Considering each Li as a linear 
map on the system of master functions Gale dual to the Laurent system F{x) = b is 

n-\-t 

___ (A\ 

(1.11) IpMs/) + W = i j = 

i =1 

When 7LA has odd index in 7L n , the solutions to F(x) = b in (R x ) n are scheme-theoretically 
isomorphic to the solutions to (11.111) in the complement of the hyperplane arrangement 
defined by Lfy) = —A, for i — 1 ,..., n+£. For any A with ZA of full rank, the solutions 
to F(x) = b in (R>o) n correspond to solutions to (11.lip where Lfy) > —A for i = 
1,..., n+A These functions (11.111) form a regular sequence in the hyperplane complement. 

The effect of the choices of B , L , and A on the Gale system and its solutions may be 
deduced from Proposition 11.11 

Corollary 1.2. Changing the basis B of the integer nullspace of A will change the Gale 
system, but not its solutions in or the hyperplane arrangement. Changing the basis L 
and the point A has the effect of an affine-linear change of coordinates in 

Remark 1.3. There are two other forms for the Gale system. First, we may take loga¬ 
rithms of the equations (11.lip to obtain the logarithmic form, 

n+t 

(1.12) := J^log(Lj(j/) +Ai) = 0 j = 1,...,£. 

2=1 

This only makes sense in the positive chamber of the hyperplane complement, 

A := {y G | Li(y) > —X, for all 7 = 1,..., n+£} . 
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To obtain the other form, set a + := max{a, 0} and a~ := max{— a, 0}, for a£Z. Note 
that a = a + — a~. With this notation, we may clear the denominators of (II.lip to get a 
polynomial form of the Gale dual system, 

n+£ n-\-£ 

(M3) n< L ‘(»)+ a <) w “ )+ - n( L ‘(»)+ A <) <e,li,) " = ° j= i, 

i =1 i =1 


Remark 1.4. It is no loss of generality to assume that this positive chamber A is bounded. 
Indeed, suppose that A is unbounded. This occurs if and only if the coefficient matrix C 
annihilates a positive vector. Since Li,, L n+ e span the dual space to A is strictly 
convex and therefore has a bounded face, P. Then there is an affine form, ao + L(y) 
where L is a linear function on and a o > 0, that is strictly positive on A and such 
that P is the set of points of A where this polynomial achieves its minimum value on A. 
Dividing by ao, we may assume that the constant term is 1. 

Consider the projective coordinate change y i-a- y, where 


(1.14) 
We have 


Vj '■= Vr 


1 + L(y) 


for j = I,...,- 


Vi = Vr 


i - m 


for j = 1, • • • 


Note that (1 + L(y))(l — L(y)) = 1. 

The coordinate change (II. 14[) manifests itself on affine forms as follows, 
affine form, then 


If p(y) is a 


(1.15) 


p(y) = p(y ) 


i - m 


and p(y) = p(y) 


i + m ’ 


where p(y) := p(y ) —p(0)L(y). If A i(y) := L^y ) + A* is one of the affine forms defining A, 
then Ai(y) = Li(y) + Aj(l — L(y)). Under the projective transformation (11.14ft , the 
polyhedron A is transformed into A, where 


A := {y G | 1 — L(y) > 0 and A i(y) >0 for i — 1,..., n+£} . 
Set A 0 (y) := 1 - L(y) and := - ES P?• 


Proposition 1.5 ([S], Prop. 2.3). Under the coordinate change (II. 14)1 . the system of 
master functions (II. lip on A is transformed into the system 


£+n 

____ 

= 1 for j = 

i =0 


on the bounded polyhedron A. 


This coordinate transformation y > y is a concrete and explicit projective transforma¬ 
tion sending the unbounded polyhedron A to a bounded polyhedron A. We will henceforth 
assume that A is bounded. 
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1.3. Khovanskii-Rolle continuation. Khovanskii-Rolle continuation is a numerical al¬ 
gorithm introduced in [5] to compute the solutions to a Gale system in the positive 
chamber. It is based on the Khovanskii-Rolle Theorem and the proof of the fewnomial 
bound [7] which bounds the number of solutions to a Gale system. 

We first state the Khovanskii-Rolle Theorem. Write Va(/i, • • •, fk) for the common 
zeroes in A of functions fi,..., fk defined on A, counted with multiplicity and J(fi,..., fe) 
for the Jacobian (determinant), det (dfi/dyf), of fa,.. ., ft. For a curve 7 in A, let ubcA( 7 ) 
be the number of unbounded components of 7 in A (those that meet its boimdary). This 
is one-half the number of points of intersection of 7 with the boundary of A. 

Theorem 1.6 (Khovanskii-Rolle |15j). Let gi,...,gg be smooth functions defined on a 
polyhedral domain A cM^ with finitely many common zeroes and suppose 7 := Va(< 7 i, •.., ge¬ 
ts a smooth curve. Let J be the Jacobian of gi,... ,gg. Then 

( 1 - 16 ) \V A (gi, ■ ■ ■, ge)\ < ubc A ( 7 ) + |V A (si, • • ■ ,gt-i, J)\ ■ 

The main idea is that along an arc of the curve 7 , the Jacobian vanishes at least once 
between any two consecutive zeroes of gg. 



If we let S be the points of 7 where J = 0 and T be the points where 7 meets the 
boundary of A, then every zero of gg on 7 lies in a unique interval of 7 \ (S U T). Thus, 
the Khovanskii-Rolle Theorem leads to the following continuation algorithm. 

Algorithm 1.7. Using arclength continuation to follow 7 starting at points of S U T, 
where 7 is traced in both directions from points of S and is traced into the interior of the 
polyhedron A from points of T, will recover all points of 7 where gg = 0, each at least 
twice. 

We apply this to a Gale system. Let fa ,..., <fg be the functions in the logarithmic form 
of a Gale system ( 11 . 121 ) . Set 'fg := V((f) 1 ,..., fa-i). Then 7 ^ is an algebraic curve in A, as 
logarithmic functions fa define the same sets as the polynomial functions (jl. 13[) . and these 
form a complete intersection. Let Jg be the Jacobian determinant of fa,...,<fg. Given 
the points Sg-\ := V(fa,... ,fa-i, Jg) where Jg vanishes on 7 ^ and the points Tg where 
7 ^ meets the boundary of A, we may compute the solutions Sg := Va(0i, • • • ,fa) using 
arclength continuation along 7 g, by Algorithm 11.71 

A genericity assumption on the coefficient matrix C and the exponents B implies that Tg 
is a subset of the vertices of A. By Lemma 3.4 (1) of [7], Jg := Jg-Ylii^iv) + Aj) is a poly¬ 
nomial of degree n. Thus we may compute the solutions Sg to the Gale system G , given 
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knowledge of the set Tj and the set 5j_i in which the transcendental function fif defin¬ 
ing Sf is replaced by the low degree polynomial J The Khovanskii-Rolle algorithm [5J 
and fewnomial bound [iTj iterate this procedure. Before explaining them, we discuss our 
genericity assumptions. 

The exponent vectors B = {/jW,... ,/3^} are the columns of a (n+£+1) x i- matrix, 
also written B , whose rows correspond to the affine functions defining A. For a face P of 
A and any j = 1 ,let B Pj be the sub matrix of B consisting of the first j columns 
of the rows of B corresponding to the affine forms vanishing on P. We assumed A is 
bounded, which may be achieved by a projective change of coordinates, as explained in 
Remark 11.41 

Assumption 1.8. We will make the following assumptions. 

(1) The affine forms Lfiy) + A* for i = 1,... ,n- fT+1 are in general position in that 
exactly i of them vanish at each vertex of A. 

(2) For every face P of A, the square matrix <Bp,codim(P) annihilates no nonzero non¬ 
negative row vector. 

(3) The sets 7 j defined below are curves. 

Assumption ( 1 ) is implied by the assumption that the coefficients of the system F(x) = b 
are general, which is mild, as we already assumed that the system had only isolated zeroes. 
A consequence is that A is a simple polyhedron. That is, any face P of codimension r 
lies on exactly r facets, and thus exactly r affine forms Lfiy) + A* vanish on P. 

Assumptions ( 2 ) and (3) depend upon the exponent vectors A of the monomials in F, 
as well as the polyhedron A (which comes from C and b ), and therefore it is more difficult 
for them to hold. As explained in [7], these assumptions hold on a Zariski open subset of 
the spaces of matrices C and B (when B is allowed to have real entries). 

We sketch the steps of Khovanskii-Rolle continuation under these assumptions. For 
each j = £, £— 1 ,..., 1 , set Jj to be the Jacobian of ..., fij, Jj+i, ■ ■ ., Ji and define 7 j 
to be Va(0i, • • •, 4>-j-\ , Jj+ 1 ,.. •, Je). By Assumption II .81 3). this is a curve. Let Tj- be the 
points where 7 j meets dA. Finally, for j = 0 ,..., £ set Sj := Va(</>i, • • •, fij, Jj+ 1 , ■ • •, Jf). 
Then Sf is the set of solutions to the Gale system. 

By Algorithm 1 1.71 arclength continuation along 7 j beginning from points of Tj and Sj~ 1 
will compute all points of Sj. Thus, the solutions Sf of the Gale system in A may be 
computed iteratively from S 0 and the sets Tj,..., Tf. Computing these sets is feasible, for 
by Lemma 3.4 of [7J and Corollary 12.31 we have 
~ / \ 

( 1 ) Ji = Ji ■ [Y[i(Li(y) + Aj is a polynomial of degree 2 e P ■ n. 

( 2 ) The points of Tj all he on some faces of A of dimension £—j and are the points of 

those faces where the polynomials Jf vanish. 

For j = 2 ,..., £, let pj := Va(4>i, ■ ■ ■, <fij- 1 )- The faces of A in ( 2 ) lie in the closure 
of the set y.j and are determined by B (Corollary 12.31 gives the precise statement). This 
set pj is algebraic, as the transcendental functions fi define the same set in A as do the 
polynomials fll,13p . Furthermore, pj has dimension £—j+ 1. Indeed, the rational function 
versions ( 11. 1 If) of the <fn{y) form a regular sequence in the hyperplane complement, so 
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their common zero set fij has expected dimension £—(j— 1). The assumption that 7 j is a 
curve is that the Jacobians Jj+i, ■ ■ ■, Je form a complete intersection on fij. 


1.4. Our running example. Let us consider this for the example of Subsection fl.il The 
first two assumptions hold. All polygons are simple, and no entry and no minor vanishes 
in the matrix of exponents, 


B = 


-1 -3 1-1 3 

2-13-2 -2 


The Jacobian of the logarithmic system (11.51) is the rational function, 

j 2x 3 — 16 x 2 y + 12 xy 2 + 6 y 3 — + 26 xy — ^-y 2 + y — 2 

xy( 1 - x - y){\-x+y)( 1 + x - 3y) 

whose denominator is the product of the linear factors defining the lines in Figure [Q 
Clearing the denominator and multiplying by 2 gives the cubic polynomial, 

J 2 := 4x 3 — 32 x 2 y + 24 y 2 x + 12 y 3 — 31x 2 + 52 xy — 53 y 2 + 9x + 15 y — 4 . 


The Jacobian J\ of g and J 2 has denominator the product of the squares of the forms 
defining the lines, and its numerator is the sextic 


- 56x 6 + 464 x 5 y - 456 x 4 y 2 - 1792 x 3 y 3 + 1896 x 2 y 4 + 336 xy 5 - 72 y 6 + 640x 5 - 1960 x 4 y 

+ 456xV + 3096 x 2 y 3 - 5176 xy 4 + 480 y 5 - 482x 4 + 2248 x 3 y - 2416 x 2 y 2 + 4256on/ 3 + 250 y 4 

+ Qlx 3 — 401 x 2 y — 1769 xy 2 — 491 y 3 — 10x 2 + 664 xy + 278y 2 — 81x — 101 y + 16 . 

The system A = J 2 = 0 has 18 solutions with 10 real, but it has only two solutions So bi 
the pentagon. The curve 71 is defined by the vanishing of the Jacobian J 2 . The set T\ of 
points where it meets the boundary of the pentagon consists of two points. 

In Figure [2] we illustrate the configuration of the sets So, SJ, S' 2 , Tf, and T 2 in the 
pentagon of Figured] along with the curves defined by Ji, J 2 (written 7 i), g (y 2 ), and /. 
The middle picture shows that arclength continuation along 7 x in both directions from 



FIGURE 2. Solving the Gale system with Khovanskii-Rolle continuation 

each point of So and into the pentagon from each point of 7\ will find the point Si twice 
(the other four continuations terminate when they either reach a point of So or exit the 
pentagon). The picture on the right shows that arclength continuation along y 2 in both 
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directions from the point of S± and into the pentagon from the four points of T 2 will reach 
the three points of S 2l twice each. These points are 

(1.17) (0.29557,0.32316), (0.14846,0.26681), (0.53346,0.20563). 

Remark 1.9. In Figured we see that some of the points T 2 are singular points of / = 0, 
but not of g — 0. In general the points T\ are smooth points on the curve 71 , but the 
other points Tj with j > 1 may be singular points of the curves 7 j they lie upon. I 11 
Section [2l we will discuss when it is possible to eliminate tracking from these possibly 
singular points. A future publication [2] detailing methods (and software) will explore 
strategies for arclength continuation along 7 from these singular points of the boundary. 

Remark 1.10. Having computed approximations to the points (11.171) . we do not yet have 
solutions to the original polynomial system. By Gale duality, if Sq is the set of solutions 
to the Gale system, then Lpf^fif^Sc)) is the set of solutions to the original polynomial 
system. Suppose Sq is a set of approximations to solutions of the Gale system. The 
difficulty is that points of iJj(Sq) likely do not lie in ^((KC)" - ), and so (fi (Sq)) is 
undefined. In Subsection 13.21 we give a method to overcome this obstruction. 


2. A Descartes’ bound for fewnomials 

We give a refinement to the upper bound of ^4^2 ( 2 ) 7 / for the number of solutions to a 
Gale system that was discovered while developing our software. This refinement involves 
sign conditions and gives a Descartes’-like bound for fewnomials. We then discuss how 
this refinement affects the running of the Khovanskii-Rolle algorithm. 

Theorem 2.1. Suppose that <3>(o:) = 0 is a system of n Laurent polynomials in n variables 
involving n +£+1 monomials. If the exponent vector of one monomial lies in the relative 
interior of the convex hull of the exponent vectors of the other monomials and if the 
coefficient matrix of these other monomials does not have a positive vector in its kernel, 
then this system has at most 


(l + 2~ l (l + ^2 (*V 


nondegenerate positive solutions. 

We will give a proof of Theorem 12.11 as part of our discussion on our improvement to 
the fewnomial bound. We first derive the fewnomial bound originally given in [7j. 

As explained at the end of Subsection II.31 the solutions Si to the Gale system are found 
by first computing So := Va(</i, ..., Jf) and certain sets Ti,..., T) where Tj consists of 
the points of the boundary of A that meet the curve 7 j. By Corollary 12.31 below, Tj 
consists of the points of certain faces of A of dimension £—j where J ]+ \,..., .If vanish. 
As Ji is a polynomial of degree 2 we may replace faces of A by their affine span and 
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use Bezout’s Theorem to estimate the cardinality of these sets, 


( 2 . 1 ) 


t 

I So | < JJdeg(Ji) = 2(2)?/ 

2—1 

£ 

\Tj\ < f,-j( A) • [] deg(Ji) = A)2fi') n '-1 , 

i=j+l 


where fi(A) is the number of z-dimensional faces of A. In fact, there is a tighter estimate 
for \Tj\ given below (12.4|) . 

From the Khovanskii-Rolle Theorem, we have \Si\ < |Si_i| + ||Tj|, and so 


( 2 . 2 ) 


N < |s„| + i(|Ti| h— + |r,|) 

< 2 «)n ! + i^/,_,(A)2(T)„<-l 



The estimate of the last row is Lemma 3.5 in [7]. 

Any difference in the two sides of any inequality (its slack ) propagates through to the 
final bound. One clear way that the estimates (12. ip could have slack would be if a system 
Jj. |_i,..., Jg on the complex affine span of a face did not have all of its solutions lying in 
that face. In practice, we have found that very few of the solutions lie on the desired face 
of A. While this seems impossible to predict or control, the Khovanskii-Rolle continuation 
algorithm always takes advantage of this source of slack. 

Another source of slack is our estimates for \Tj\, which depend on the number of certain 
faces of A, as not all faces of codimension j can contain points of Tj. We investigate the 
location of the points Tj. For each j — 2, we defined jij C A to be closure of the 

points in the interior of A given by 


(2.3) 0i (y) = ••• = (pj-i(y) = 0. 

This set is algebraic and has dimension £—j+ 1. 

Lemma 2.2. Suppose that Assumption 11.81 (1) holds. If P is a face of A whose interior 
meets pj, then the submatrix Bp,j-i of B consisting of the first j —1 columns of rows 
corresponding to the facets of A containing P annihilates a nonzero nonnegative row 
vector. When this occurs, we have P C /i,. 


By annihilating a nonzero nonegative row vector, we mean there is a nonzero vector 
v = (wi,..., v r ) with Vi > 0 for all i such that vBpj_ 1 = 0. Here, r is the number of facets 
containing P. 
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Proof. Let p be a point lying in the interior of P. After an affine change of coordinates 
and a reordering, we may assume that p — 0 and the affine forms defining A are 


yi, V2, ■■■, y r , K+i(y) , • • •, A n+t+iiy), 

where A a (0) > 0 for a > r. That is, points of P have their first r coordinates equal to 0. 
Then pj is defined on the interior of A by the rational functions 


(0 


Pi 

Vl ■ ■ ■ y'r 


q 0) 


n+i+1 


A a(y)^ = 1 for 2 = l,...,j-l. 


a=r-\-1 


After an analytic change of coordinates in the neighborhood of the origin p, we may 
assume that these become monomial, 

o(i) 

yi ••• yr r = 1 for i = 1, - - -, j—1 - 


The solution set meets the origin p if and only if it contains a monomial curve of the form 


(■ t ai , t 


Ot-2 


, t ar ) 


with each exponent cq nonnegative and not all are zero. But then a = (oq,..., ot r ) satisfies 

r 

^= 0 for i = 1 ,..., j-1 , 
t= i 


so that a is a nonzero nonnegative row vector annihilated by Bp,j- 1 - 

The last statement follows as these arguments may be reversed. □ 


Corollary 2.3. Suppose that Assumption 11.81 holds. Then pj fl dA is the union of those 
codimension j faces P of A where Bp.j-\ annihilates a nonzero nonnegative row vector. 

The points Tj are the points of pj D dA where the polynomials J J+ \,..., J( vanish. By 
Corollary 12.31 this is the restriction of those i—j polynomials to a collection of faces, each 
of dimension i—j. 


Proof. Since pj has dimension i—j A- 1, its intersection with the boundary of A has dimen¬ 
sion i—j. If this intersection meets the relative interior of a face P of A of codimension 
r, then the matrix Bpj -1 annihilates a nonzero nonnegative row vector, by Lemma [2.21 
If r < j, then Bp, r is a submatrix of Bpj- 1 , and so annihilates a nonzero nonnegative 
row vector and contradicts our assumption. Thus r > j, and we see that the rest of the 
corollary follows from Lemma 12.21 □ 

Corollary 12.31 restricts the faces on which points of Tj can lie and implies the following 
tighter estimate, 

(2.4) \Tj\ < M e _j(A,B) -2 e - j n £ - j for j> 1, 

where Mi_j(A, B) counts the faces P of A with codimension j for which annihilates 

a nonzero nonnegative row vector. 
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Proof of Theorem \2.1[ Let $(x) = 0 be a system of n Laurent polynomials in n variables 
involving n+£+l monomials. Let x° be a monomial in <f> whose exponent vector a lies 
in the interior of the convex hull of the other exponent vectors. We may replace <f> by 
x~ a Q to obtain a system of Laurent polynomials with the same zeroes in (M x ) n having a 
constant term such that 0 lies in the interior of the convex hull of the exponent vectors 
of its remaining monomials. 

Rewrite $(x) =0 as F(x) = b where F has n+£ monomials with exponent vectors A, 
and the origin lies in the interior of the convex hull of A. Thus the origin is a nonnegative 
integer combination of the vectors in A, with every vector in A participating in this 
expression. Consequently, this expression is a positive vector in the integer nullspace of 
A , and we may assume it is the first column of the matrix B of exponents for the Gale 
system. This implies that no j x (j—1) submatrix of the first j —1 columns of B has a 
positive row vector in its kernel, and so M^_j( A, B) — 0, when j > 1 . 

If the positive chamber A is bounded, then it has at most n+£ facets and T 2 ,..., T) are 
empty. This gives the estimate 

\S t \ < |-So| + ||Ti| < 2 ( 2 ) 7 /+(n+£)2(^“V- 1 

(2.5) = (l + 2~ l (l + ^2 ( 2 ) 7 /. 

If A is unbounded, then, after the transformation of Remark 11.41 A possibly has n+Z+1 
facets and the first column of B has exactly one negative entry. 

To complete the proof, suppose further that the coefficient matrix C of those other 
exponent vectors has no positive vector in its kernel. This is equivalent to A being 
bounded, and we have the bound (12.5[) . □ 

The sets T 2 ,..., T) are empty under the hypotheses of Theorem 12.11 and when the first 
column of B has every entry positive. In this case, Assumption 11.81 ( 2 ) on B is satisfied. 

We noted in Remark II.91 that the points Tj for j > 2 will typically be singular points of 
the curve 7 j as ji 3 is singular along the boundary of A. Minimizing the number of points Tj 
(or eliminating them altogether as in Theorem 12.ip . reduces the difficult path tracking 
required to track from these singular points Tj, thereby improving the performance of the 
Khovanskii-Rolle continuation algorithm. 

3. Algorithmic details for the Gale transform 

Obtaining the solutions of a fewnomial system by Khovanskii-Rolle continuation is a 
3-step process. 

(1) Apply the Gale transform to a fewnomial system F(x) = b, producing a dual Gale 
system G(y) = 1 of master functions. 

(2) Use Khovanskii-Rolle continuation to find approximations to elements of Sq which 
are the real solutions of G(y) = 1 in the positive chamber Acl f . These approx¬ 
imations form the set Sq. 

(3) For each s* G Sq, transform s* to t*, a numerical approximation to t G Sf- Set 
Sp := {t* | t* is an approximation to t G Sf}. 
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A fewnomial system may be solved using Gale duality by other root-finding methods 
in place of Khovanskii-Rolle continuation in Step 2. We provide pseudocode for the 
Gale transform of Step 1 (Algorithm 13.1ft and for the transformation of the approximate 
solutions Sq of the Gale system G to the approximate solutions Sp of the fewnomial 
system F in Step 3 (Algorithm 13.3p . Two steps of Algorithm 13.11 require extra attention; 
they are detailed in the final two subsections. 

3.1. Converting a fewnomial system to a system of master functions. Using the 
notation of Section [fl we give pseudocode for the basic algorithm for converting a system 
F(x) = b of fewnomials into a system of master functions G(y) = 1. 

Algorithm 3.1. 

Input: £, n, an n x (n+£) integer matrix A , an n x (n+£) coefficient matrix C, and target 
column vector b. 

Output: A system G(y) = 1 of l master functions in £ variables that is Gale dual to the 
system F(x) = b. 

1) Choose £ monomials and reorder the columns of C so that these monomials are the 

last £ columns. 

2) Put [C | b] in echelon form, obtaining a matrix [— I n \ L | A]. 

3) From the rows of L and entries of A form the affine forms L^y) + A,. 

4) Choose l integral basis vectors /% of the nullspace of A to construct B = {^\ ■ ■ ■, /3^}. 

5) Use B and the forms Li(y) + A i to form the Gale dual system G(y) = 1 as in (11,lip . 

We explain some of details in the design and implementation of each step of Algo¬ 
rithm 13.11 and then illustrate it with an example. 

1) The effect of this choice is an affine change of coordinates in R , which affects the 

shape of the positive chamber and the performance of the Khovanskii-Rolle 

continuation when no scaling of the Gale system functions occur. With a simple 
scaling routine, discussed in Section [4] the parameterization of the monomials 
appears to not matter when £ = 2. 

2) We simply use Gaussian elimination with partial pivoting. 

3) If the space has coordinates xel" and i/Gl*, the echelon form [— I n \ L | A] 
gives a parameterization of the set C(z,y ) = b by the affine forms Zi = L t (y) + A* 
for i — 1 ,..., n. 

4) We use the LLL algorithm [17] to compute an initial basis B and then find vectors in 
the row space with few negative entries, a heuristic to minimize A, B) (12.4j) . 
which is the number of codimension j faces of A on which the points Tj may lie. 
This step is the focus of Section 13.31 

5) If we set L n+ j(y) := yj and A n+ j := 0 for j = 1,...,£, then the vectors in B 
and affine functions Li(y ) + Aj for i — 1 ,... ,n+£ together give the Gale system 
G(y) = 1 as in (II. lip in Proposition 11.11 
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Example 3.2. Consider the fewnomial system, where n = 5 and 1 = 2 



-a 1 b 2 c 2 d 

+ 

{\b 2 c 

+ 

2 b~ 4 c~ 

d~ 5 e~ 

-1 

— 

1) 

= 0 , 


—ac 

+ 

\{\b 2 c 

— 

b~ 4 c 

~ 7 d~ 5 e~ 

-1 

+ 

1) 

= 0 , 

(3.1) 

— be 4 d 4 

+ 

\{~\b 2 c 

— 

3 b~ 4 c 

~ 7 d~ 5 e~ 

-1 

+ 

6 ) 

= 0, 


— d 

+ 

\{~\b 2 c 

— 

2 b~ 4 c 

~ 7 d~ 5 e~ 

-1 

+ 

8 ) 

= 0, 


—e 

+ 

{-\b 2 c 

+ 

2 b~ 4 c 

~ 7 d~ 5 e~ 

-1 

+ 

3) 

= 0, 

with support {a~ 

l b 2 c 2 d , ac , bc 4 d 4 , d , 

e, l 

9 c , lr 

4 c~ 7 d- 

5 e- 

- 1 } 

and thus matrix of expo- 


nents 


and coefficient matrix 


[c I b] = 



-1 

1 

0 

0 

0 

0 

0 


2 

0 

1 

0 

0 

2 

-4 

A = 

2 

1 

4 

0 

0 

1 

-7 


1 

0 

4 

1 

0 

0 

-5 


0 

0 

0 

0 

1 

0 

-1 


' -1 

0 

0 

0 

0 

1/2 

2 

-1 ' 

0 

-1 

0 

0 

0 

1/8 

-1/2 

1/2 

0 

0 

-1 

0 

0 

-1/16 

-3/4 

3/2 

0 

0 

0 

-1 

0 

-3/8 

-1 

4 

0 

0 

0 

0 

-1 

-1/2 

2 

3 


where the first seven columns correspond to the seven monomials in A and the final 
column corresponds to the constants. 

We presented this fewnomial system with a coefficient matrix in echelon form for the 
given order of the monomials, so Steps 1 and 2 of Algorithm 13.11 are complete. If we 
set s := b 2 c and t := b~ 4 c^ 7 d~ b e ~ 1 , the rows of this echelon matrix express each of the 
first five monomials as affine functions of the parameters s and t, so that, for example, 
a~ l b 2 c 2 d = + 2t — 1. We have the seven (=5 + 2) affine forms, 


(3.2) 


Ai(s, t) = |s + 2t — 1 
A 2 (s,t) = |s — \t + \ 

A s(M) = -^-|f+| 
A 4 (s,f) = -§s-t + 4 
A 5 (s, t) = — 7 + + 2t + 3 
A 6 (s,t) = s 
A 7 (s,t) = t 


Given a vector /? = [/fi,..., /3 7 ] T in the nullspace of A, we have ]/[J = i(x ai f * = where 
di is the i th exponent vector in A. Under the substitution x ai = A i(s,t), we obtain 
f | i=1 Aj(s,f)A = 1. One choice of basis B for the nullspace of A is 

{[4,4, 2 , 3, 3,1, 3] t , [-1 ,-1, 2,-2,1,2, If}. 
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This gives the system of master functions (written in the polynomial form (jl.l3]i ) 

g(s,t) = (s+4f—2) 4 (s—4f+4) 4 (s+12f—24) 2 (3s+8f—32) 3 (s—4f—6) 3 st 3 
(3.3) -68719476736, 


f(s,t) = —2(s+12t—24) 2 (s—4t—6)s 2 f — (s+4t—2)(s—4t+4)(3s+8f—32) 2 . 

Figure [3] displays the algebraic curves in R 2 defined by these polynomials, together with 
the lines defined by the affine forms Aj(s, t) = 0 for % = 1 ,..., 7. 



Figure 3. Master function curves in the septagon 

The original system (13.ip has 95 complex solutions with 11 real and six in the positive 
orthant. The hypotheses of Theorem 12.11 hold: with n = 5 and i = 2, its improved bound 
is 87 by Theorem 12.11 and 67 using (12.2ft . These hypotheses imply that the curve g — 0 
does not meet the boundary of the septagon A. 

3.2. Transforming Gale system solutions to fewnomial solutions. The output of 
Khovanskii-Rolle continuation is a set S* G of approximations to solutions S G of the Gale 
system in the positive chamber A. Converting these points of A to points Sp of R> that 
are approximations to solutions Sp of the fewnomial system is problematic, for there is 
no natural map between points of A and R”, except on the solutions to the two systems. 
We present a simple method to accomplish this conversion. 

Let s* G Sq and suppose that s is the point of Sq that s* approximates. Applying 
the map -0 (11 ,8ft of Subsection 11.21 whose components are the affine functions Li(y) + A* 
gives points 0(s*) and 'ip(s) in R> +£ that lie on the affine plane H . The point 0(s) lies on 
^(M>), but the point ip(s*) does not lie on ^(R”) (unless s* G Sq), as the former is a 
solution of the Gale system, while the latter typically is not. 

To remedy this, choose a square submatrix A! of rank n of the matrix A of expo¬ 
nents. Projecting M. n+e to those coordinates indexed by the columns of A 1 gives a map 
7 ^ 4 /: M n+f —>• M" . The composition 

R” -HA-} AAA. 
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is the map 4 /, which is invertible on R™. We set 

t* ■= o n A , o 7p(s*). 

Since t = pT^, o 774 / o -0(s) = p^ o 0(s) is the solution to the fewnomial system, and these 
maps are differentiable, we expect that t* is close to t. We formalize this transformation. 

Algorithm 3.3. 

Input: The exponent matrix A and affine functions Li(y) + A* from Algorithm 13.11 and 
the set Sq = {s* | % — 1 ,..., A} of approximations to solutions of the Gale system G in 
the positive chamber A. 

Output: A set Sp = {t* \ i = 1,...,1V} of approximations to positive solutions of the 
original fewnomial system F. 

Choose a full rank n x n sub matrix A 1 of the matrix A and compute (A') -1 , 
for % := 1 to N do 

a) Evaluate the map 0 at the point s* to get the point z* G R> + ^, 
z i := (-^1 ( s i) + ■ ■ • j L n+e (s*) + \ n+ i) . 


b) Project z* to R> using set v* := n a'{z*). 

c) Let w* e M n be the coordinatewise logarithm of v*. 

d) Set x* := {w^A 1 )" 1 ) 7 ^ and let t* be the coordinatewise exponentiation of x*. 

end for 

Note that the homeomorphism ip. 4 /: R> —> R” becomes the linear isomorphism induced 
by x i-A- A'x in logarithmic coordinates. Steps c) and d) invert this transformation. 

Even if s* is certified to be an approximate solution to the Gale system with associated 
solution Si (i.e. s* —$■ s* under Newton iterations, doubling the number of significant digits 
with each iteration), then it need not be the case that t* is an approximate solution 
to the fewnomial system with associated solution L := (p~^ o -0(sj). However, we may 
check this using, for example, the software alphaCertified |13j . If t* is not certified 
to be an approximate solution, then we could refine s* to be closer to s* and reapply 
the transformation. Our software has the functionality to do this, either producing a 
soft certificate (all calculations done in floating point arithmetic) or a hard certificate 
(calculations done with rational arithmetic), which is a proof that t* is an approximate 
solution. 

Example 3.4. An approximation to a solution of f(s,t ) = g(s,t) — 0 from Example 13.21 
lying in the positive chamber is s* = (0.94884808,0.65721633). We follow the steps of 
Algorithm 13.31 to recover a solution of the original fewnomial system. 
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We first choose an invertible sub matrix of A. The first 5(= n) columns of A suffice as 
these are linearly independent. Evaluating the map i/j at gives the point 

' 0.78885671 ' 

0.28999784 

0.94778474 

z* = 2.98696564 . 

3.84000863 

0.94884808 

0.65721633 


Projecting to the first 5 coordinates of z* and taking coordinatewise logarithms yields 


-0.23717058 

-1.23788180 

-0.05362787 

1.09425803 

1.34547461 


Proceeding to the next step, we evaluate 


x 


* 

1 


{w* 1 {A!~ 1 )) T 


0.020520123 

0.60294767 

-1.25840192 

1.09425803 

1.34547461 


Finally, 


t\ = exp(a^) 


1.0207321 

1.8274977 

0.28410769 

2.9869656 

3.8400086 


We carry out Algorithm 13.31 with sj' starting with 135 bits of precision for the input and 
truncate the output to an approximation that can be stored in one 64-bit floating point 
number. 

Using alphaCertif ied [13], we perform two Newton iterations on this approximation, 
doubling the precision each time. We obtain a soft-certification of this refined approxi¬ 
mation, and then convert the numerical approximation to its rational form and obtain a 
hard certification that proves that t* converges to t\. 


3.3. Heuristic for the nullspace basis. A choice in the Gale transform that affects the 
efficiency of numerical tracking in the Khovanskii-Rolle algorithm is the selection of the 
basis B of the annihilator of A. 

Observe that the points of Tj for j > 2 are often singular points of the curve 7 j. For 
efficiency and numerical stability, we want to avoid tracking near these points. Corol¬ 
lary 12.31 implies that the curve 7 j may approach a codimension j face P of A only when 
Bpj -1 annihilates a nonzero nonnegative row vector. By Corollary 11.21 changing the basis 
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elements /3^ of B will change the Gale system as well as the curves 7 j, without changing 
the solutions. Ideally, we would use this freedom to choose the basis B of the annihilator 
of A to minimize the number of faces of P on which points Tj for j > 2 could he. 

The number of possible choices for B up to patterns for the signs of the coordinates 
suffers from a combinatorial explosion [T4j as i increases. Consequently, this optimization 
problem may become infeasible for i large. Instead, we find elements of the annihilator 
of A that have the fewest number of negative coordinates. We expect that this choice 
will tend to minimize the number of faces on which the points Tj could lie. In the case 
of Theorem 12.11 this choice achieves the minimization. When there is a positive vector 
annihilating A, we may choose B to consist only of positive vectors, and so no submatrix 
annihilates any nonzero nonnegative row vector. By Lemma [2721 no face of P will contain 
points of Tj. 

Once a basis for the annihiliator of A is known, the algorithm proceeds differently 
depending upon the size of l. For small £, we can enumerate all possible patterns for 
the signs of the coordinates of vectors in the annihilitor of A, and then choose a basis B 
consistsing of vectors having the fewest number of negative coordinates. When i is large, 
we take random integer linear combinations of a basis of the annihilator of A and choose 
of these generated vectors those with the minimum number of negative entries. 

This heuristic does not always work in parctice, as we observed in our test suite de¬ 
scribed below. One problem is when a point of some Sj —a solution to an intermedi¬ 
ate problem—lies close to the boundary of A. Such solutions may be numerically ill- 
conditioned, and the subroutine in the Gale duality software package that tracks the 
curves 7 , near the boundary of A (the monomial tracker described in [5], § 4.3]) may 
miss such ill-conditioned solutions. A future publication [2] will focus on the design and 
performance of the monomial tracker subroutine. 

4. Software 

The algorithms of Section [3] are implemented in a software package, galeDuality [3], 
which is available at each author’s website and written in C++. This will be the front end 
for a package we are developing for Khovanskii-Rolle continuation. We describe a test 
suite (14.11) to evaluate the proposed heuristic in Section 13.31 and some implementation 
details of the package (14.21) . 

4.1. Test suite. As noted in Sections [T] and [2 Khovanskii-Rolle continuation finds the 
solutions Si to a Gale system in a polyhedron A using arclength continuation along 
curves 7 j starting at points T 0 := S 0 , 7),..., T f where Tj are solutions to £ — j Jacobian 

determinants (having Bezout number 2^ 2 in certain codimension j faces of A. The 
set To lies in the interior of the polyhedron A and points of T\ could lie in every facet, 
but the other Tj lie in only a subset of the faces of codimension j. The heuristic of 
Subsection 13.31 is intended to minimize or eliminate the number of faces contributing to 
Tj for j > 2 . 

To test this heuristic, we created a test suite of Gale systems for each (£, n) pair with 
i = 2 and n — 2,..., 13. 
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For each n, we generated an (n+2)-gon by successively selecting random rational points 
in an annulus until we had n +2 points in convex position. These points formed the vertices 
of the (n+ 2 )-gon. The edges of this polygon then give n+2 affine forms with relatively 
prime integer coefficients. To create a Gale system, we generated two vectors of exponents 
whose components correspond to the affine forms. Both had integer components with 
absolute values selected uniformly at random in the range 1,..., 10. One vector, /3 + , had 
all components positive, but the other, j3 ± , had signs alternating as much as possible 
traversing the polygon cyclically. This was intended to emulate the configuration of the 
master function curves in Figure [3j 

Our test suite consists of 100 random Gale systems for each n. For each random system, 
we conducted two types of runs of Khovanskii-Rolle continuation, one for each ordering 
of the exponent vectors j3 + and /3 ± as columns in B. When (3 + is the first column of B, 
the set T 2 is empty, by Theorem 12.11 but when /J ± is the first column, T 2 consists of the 
vertices defined by affine forms such that the components of /5 ± had differing signs. Our 
proxy to test for the affect of the heuristic in Section 13.31 was to compare the average of 
10 running times for these two different choices. 

For any fixed n, in a head-to-head comparison of the average run-time of Khovanskii- 
Rolle continuation for the two different choices, choosing /3 + as the first column of B 
outperformed /3± for > 92% of the systems in the test suite. For the majority of the 
instances in which the choice of /5 ± outperformed (3 + , the difference between the run 
times was insignificant. 

In about 15% of all the systems in the test suite, the choice of /3 + in the current imple¬ 
mentation of the Khovanskii-Rolle continuation algorithm yielded incorrect results. This 
occurs because some curve 7 j is close to the boundary of A. The monomial tracker subrou¬ 
tine fails to find some point Sj when Sj is on the curve 7 j when 7 j is being approximated 
as a monomial curve near the boundary of A. As n increased, we observed that this 
event occurred more frequently. We believe this is due to the inadequacy of the monomial 
tracker in the proof-of-concept implementation released with [5], and improving that is 
a future project. 

In addition, since the curves 7 j are defined by the vanishing of a combination of Jacobian 
determinants and master functions, all these functions may have large coefficients, as 
seen in the examples we gave. This affects the efficiency of the arclength continuation 
and computation of the points Tj. Using the logarithmic form of the master functions 
mitigates their possible contribution to this inefficiency. For the Jacobian factors, we 
propose a simple scaling heuristic that is a naive version of SCLGEN [18] : divide each 
Jacobian by the average of the largest and smallest absolute values of its coefficients. 
Scaling serves a secondary role in that it is sometimes necessary to precondition the 
systems that precompute the points Tj. Without this scaling, we were not always able to 
compute some of the starting points that lie in the boundary of A in some of the examples 
of our test suite. 

Timing of the Khovanskii-Rolle continuation algorithm, information on the number of 
points found in Tj, and the number of real solutions for each of the systems in the test 
suite is archived with our software, galeDuality [3j- 
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4.2. galeDuality. Our software package, galeDuality, is an open-source program writ¬ 
ten in C++ and accepts as input either a fewnomial system or a Gale system. When given 
a fewnomial system, the software symbolically computes the Gale transform and saves the 
Gale system to a hie, using the algorithms of Section [3l for arbitrary l. Similarly, when 
given a Gale system, galeDuality will transform the system of master functions into its 
dual fewnomial system and save the fewnomial system to a hie. 

if the polyhedron A is unbounded, the software applies the projective transformation 
of Proposition 11.51 When i = 2, galeDuality calls the maple script of [5] to solve the 
Gale system using Khovanskii-Rolle continuation. 

After the Khovanskii-Rolle continuation computation is finished, galeDuality will con¬ 
vert the approximations to solutions of the Gale system back to approximations to so¬ 
lutions of the fewnomial system using Algorithm 13.31 and then call alphaCertif ied, 
if desired by the user, to guarantee that the approximations converge to solutions. If 
the user wishes, galeDuality will also certify the approximations in A found by the 
Khovanskii-Rolle continuation algorithm. 

Finally, software and a user manual with more detailed information is available online 
from each of the authors’ websites. 
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